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ABSTRACT 

X-ray emitting clusters of galaxies are considered in the context of modified New- 
tonian dynamics (MOND). I show that self-gravitating isothermal gas spheres are not 
good representations of rich clusters with respect to the radial gas density distribution 
as indicated by the X-ray surface brightness. Pure gas spheres with a density distribu- 
tion described by a model, as observed, also fail because, with MOND, these objects 
are far from isothermal and have a gas mass and luminosity much larger than observed 
for clusters of the same mean temperature. These problems may be resolved by adding 
an additional dark mass component in the central regions; a constant density sphere 
contained within two core radii and having a total mass of one to two times that in the 
gas. With this additional component, the observed luminosity-temperature relation for 
clusters of galaxies is reproduced. When the observed X-ray surface brightness distri- 
bution in actual clusters is modeled by such a two-component structure, the typical 
mass discrepancy is three to four times smaller than with Newtonian dynamics. Thus 
while MOND significantly reduces the mass of the dark component in clusters it does 
not remove it completely. I speculate on the nature of the dark component and argue 
that this is not a fundamental problem for MOND. 
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1 INTRODUCTION 

The modified Newtonian dynamics (MOND) is an 
empirically-based modification of Newtonian gravity or in- 
ertia in the limit of low accelerations (< Qq ~ cHo) sug- 
gested by Milgrom (1983) as an alternative to cosmic dark 
matter. In addition to explaining galaxy scaling relations 
(Tully-Fisher, Fab er- Jackson, Fundamental Plane), this sim- 
ple algorithm allows one to accurately predict the shapes 
of spiral galaxy rotation curves from the observed distribu- 
tion of gaseous and stellar matter. MOND also accounts 
for the kinematics of small groups of galaxies (Milgrom 
1998) and of superclusters, as exemplified by the Perseus- 
Pisces filament (Milgrom 1997) without the need for un- 
seen mass. These well-documented phenomenological suc- 
cesses (Sanders & McGaugh 2002 and references therein) 
challenge the cold dark matter (CDM) paradigm and pro- 
vide some support for the suggestion that the current theory 
of gravity and inertia (General Relativity) may need revision 
in the limit of low accelerations or field gradients. 

However, problems do arise when one attempts to ap- 
ply MOND to the large clusters of galaxies. The and White 
(1988) first noted that, to successfully account for the dis- 



crepancy between the observed mass and the traditional 
virial mass in the Coma Cluster, the MOND acceleration pa- 
rameter, supposedly a universal constant, should be about a 
factor of four larger than the value implied by galaxy rota- 
tion curves. With MOND, the dynamical mass of a pressure 
supported system at temperature T is M oc T^/flo; there- 
fore, the The and White result could also be interpreted as 
an indication that the MOND dynamical mass is still larger 
than the detectable mass in stars and gas. 

In astronomical tests involving an individual extragalac- 
tic object, such as the Coma cluster, a contradiction is not 
necessarily a falsification. One can always argue that the 
peculiar aspects of an object, such as deviations from spher- 
ical symmetry or incomplete dynamical relaxation, exempt 
that particular case. However, Gerbal et al. (1992), looking 
at a sample of eight X-ray emitting clusters, noted that the 
problem is more general: although MOND reduced the New- 
tonian discrepancy by a factor of 10, there is still a need for 
dark matter, particularly in the central regions. Later, con- 
sidering a large sample of X-ray emitting clusters, I found 
that the mass predicted by MOND remains, typically, a fac- 
tor of two or three times larger than the total mass observed 
in the hot gas and in the stellar content of the galaxies 
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(Sanders 1998). More recently, Aguirre, Schaye & Quataert 
(2001) pointed out that MOND is inconsistent with the ob- 
served temperature gradient in inner regions of three clusters 
for which such data is available. Again, the problem can be 
remedied by additional non-luminous mass, primarily in the 
inner regions, of order two or three times the observed gas 
mass. This discrepancy is also evident from strong gravita^ 
tional lensing in the central regions of clusters- the format 
tion of multiple images of background galaxies. Here, MOND 
docs not apply because accelerations arc Newtonian, and the 
implied surface density greatly exceeds that of visible mat- 
ter and hot gas (Sanders 1998). So, although MOND clearly 
reduces the classical Newtonian mass discrepancy in clus- 
ters of galaxies, there still remains a missing mass problem, 
particularly in the cores. 

Here 1 consider the issue of the remaining missing mass 
in clusters and whether or not this a fundamental problem 
for MOND. First I calculate the structure and X-raj^ surface 
brightness distribution of MOND isothermal gas spheres. 
Except for the very central regions, the structure of these 
objects is self-similar; the finite mass is primarily determined 
by the temperature and is very weakly dependent upon the 
central density. Thus there exists a mass-tcmpcraturc rela- 
tion (M oc T^) which is absolute; this implies a gas mass typ- 
ically 5 to 10 times larger than that observed in X-ray emit- 
ting clusters of the same temperature. Moreover, for central 
electron densities in the range of 0.001 to 0.01, there is a well- 
defined X-ray luminosity-temperature relation which is less 
steep (1/ oc T^'^) and lies well above the observed luminosity- 
temperature relation for clusters (L oc T^'®). Looking at 
individual objects, the radial dependence of X-ray surface 
brightness does not reproduce that typically observed in X- 
ray emitting clusters- observations which are well-fit by the 
traditional "/J- model" (Sarazin 1988). The conclusion is that 
self-gravitating MOND isothermal gas spheres are not good 
representations of clusters of galaxies. 

Gas spheres with a density distribution described by a 
/3- model are not isothermal in the context of MOND. A cen- 
tral or boundary temperature must be specified to determine 
the run of temperature in such objects, but those models 
with the lowest temperature gradients have large core radii 
and are again over-massive and over-luminous with respect 
to observed clusters of the same temperature. These prob- 
lems may be solved by postulating the existence of a second 
rigid component in the mass distribution: a constant density 
core having a radius about twice that of the (3 model core 
and a central surface density comparable to ao/G. Here, by 
rigid I mean a component with a fixed density distribution 
which does not respond to the gravitational field of the hot 
gas or galaxies. The presence this additional component in 
the Coma cluster is implied by the MOND hydrostatic gas 
equation for cumulative mass. If this component is generally 
present in clusters, it contributes to the gravitational force in 
the inner regions and reduces the core radius at a given tem- 
perature. In this way the observed luminosity-temperature 
and mass-temperature relations for clusters may be repro- 
duced. 

About 40 individual clusters have been considered in 
terms of such a two-component model; i.e., the observed 
surface-brightness distributions and mean temperatures are 
fit by specifying the density of the non-luminous component 
which is assumed to extend to two gas core radii. The to- 



tal mass of this additional component varies between a few 
times 10^^ and 10^'* M© and the implied mass-to-light ratio 
is typically 50 in solar units. Therefore, the required rigid 
component is not a standard stellar population; as noted 
previously, MOND requires dark, or heretofore undetected, 
matter in the central regions of rich clusters. Although the 
total discrepancy between dark and detectable mass is re- 
duced by, on average, a factor of four over that implied by 
Newtonian dynamics, it is clear that there remains a dark 
matter problem for MOND. I discuss the issue of whether or 
not this is a contradiction. 1 speculate on the possible nature 
of this non-luminous component and argue that neutrinos of 
finite mass are a possible candidate. 



2 THE STRUCTURE AND PROPERTIES OF 
MOND GAS SPHERES 

2.1 MOND isothermal spheres 

The structure of isotropic isothermal spheres may be deter- 
mined by solving the equation of hydrostatic equilibrium: 



dp 
dr 

with the pressure p given by 



-pg 



p = par 



(la) 



(16) 



where ar is the one-dimensional velocity dispersion (con- 
stant for an isothermal sphere), p is the density. The grav- 
itational acceleration, g, in the context of MOND, is given 
by 



gp.{g/ao) = g„. 



(2a) 



where gn is the usual Newtonian gravitational acceleration. 



gn 



GM(r) 



(26) 



flo is the critical acceleration below which gravity deviates 
from Newtonian (found to be about 10~** cm/s~^ from fits 
to galaxy rotation curves), and ij,{x) is a function which 
interpolates between the Newtonian regime {p.{x) = 1 when 
a; >> 1) and the MOND regime {p{x) = x when x « 1). 
A function having this asymptotic form. 



/^(a;) = x{l + x ) 



(2c) 



works well for galaxy rotation curves and is also used here. 

Milgrom (1984) demonstrated that MOND isothermal 
spheres have a finite mass and a density which falls of as 
r~" where q ~ 4 at large radii. In the outer regions, where 
M(r) = M = constant, it follows immediately from eqs. 1 
and 2 that 



(3) 



which means that the mass is uniquely determined by the 
specified velocity dispersion or temperature. For an isother- 
mal gas sphere this relation becomes 



M 



(2.9 X lO")Tfeey'M0 (4) 



16 / kT \ 
Gtto \ frup J 

where / is the mean atomic weight (« 0.62 for an fully 
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Figure 1. The gas mass-temperature relation for clusters of 
galaxies. The points are measured temperatures and inferred gas 
masses for the 42 clusters listed in Table 1. The gas mass is that 
obtained by extrapolating the /3-model to a radius where the 
gas density falls to 10"'^* g/cm^. The dashed line is the mass- 
temperature relation for MOND isothermal spheres (eq. 4) and 
the solid curve is the gas mass-temperature relation for the two 
component MOND /3 models discussed in Section 3. 



ionized gas with solar abundances) and is the proton 
mass. 

This would be, in effect, the extension of the Faber- 
Jackson relation to clusters of gala^cies (Sanders 1994). How- 
ever, the observational definition of such a relation is am- 
biguous because the total gas mass defined by the /3-model 
is typically divergent. If one considers the mass inside a fixed 
radius (Sanders 1994), or within a radius where the density 
falls to some fixed value (Mohr, Mathiesen & Evrard 1999), 
a relationship of this form (eq. 4) is observed for clusters. 
However, the mass at a given temperature is typically a fac- 
tor of 5 to 10 below that implied by eq. 4. This is shown 
in Fig. 1 where the MOND mass-temperature relation for 
isothermal gas spheres (dashed line) is compared to the ob- 
served gas mass-temperature relation for 42 clusters listed 
in Table 1. Here the observed gas mass is given by the (i- 
model (eq. 7 below) extrapolated to a radius where the gas 
density has fallen to about 10~^* g/cm^, or 250 times the 
mean cosmological density of baryonic matter. 

The characteristic scale of MOND isothermal gas 
spheres is 

Tm « CTr^ /ao (5) 

(Sanders & McGaugh 2002). Given that the X-ray luminos- 
ity due to free- free emission is L oc rie^T^rm'^; then, eqs. 4 
and 5 above would imply an X-ray luminosity-temperature 
relation of the form L oc T The observed luminosity- 
temperature relation (e.g., Ikebe et al. 2002) is significantly 
steeper, L oc T^'^, than that of MOND isothermal spheres. 
To pursue this in more detail, I calculated the X-ray lu- 



Figure 2. The X-ray luminosity-temperature relation for clusters 
of galaxies. The points are the clusters tabulated by Ikcbc ct al. 
scaled to h=0.7. The dashed curve is the relation for MOND 
isothermal spheres and the dotted curve is the relation for the 
neaJT-isothermal MOND /3 models. The solid curve is the relation 
for the MOND two-component models discussed in Section 3. 



minosity for MOND isothermal spheres with electron densi- 
ties and temperatures similar to those of the X-ray emitting 
clusters of galaxies. I assume a central electron density of 
0.006 cm~^ as being typical of clusters, and, for a given tem- 
perature, numerically integrate eqs. 1 and 2 from the center 
outward. The optically thin free-free radiation for the entire 
sphere is then calculated from the run of electron density. 
The resulting luminosity-temperature relation is shown by 
the dashed curve in Fig. 2. compared to the observations 
by Ikebe et al. (2002) scaled to h=0.7. In both cases, this 
is the radiation emitted in the band 0.2 to 2.4 keV where 
the X-ray flux is typically measured by satellites such as 
ASCA and ROSAT. Not only is the theoretical dependence 
shallower than observed, the predicted luminosities are an 
order of magnitude larger than actual clusters at the same 
temperature. 

On the basis of these scaling relationships, it would seem 
that MOND isothermal spheres axe not good representations 
of clusters of galaxies. This conclusion is reinforced when 
we consider the surface brightness distribution of a MOND 
isothermal gas sphere. Although the mass of an isotropic, 
isothermal sphere is effectively determined by the tempera- 
ture, the detailed structure depends upon the central den- 
sity; for a single temperature there is a family of solutions 
bounded by a limiting solution with a 1/r density cusp at 
the center (Milgrom 1984). The limiting MOND solution re- 
sembles the Hernquist model (Hernquist 1990) with the 1/r 
cusp steepening to a 1/r'' beyond a radius r-m- Lower density 
spheres are characterized by a constant density core with a 
l/r'* gradient at large r; the central gas densities observed in 
clusters of galaxies would place these objects in this category 
of the sub-limiting solutions. The X-ray surface brightness 
distribution resulting from such a MOND isothermal sphere 
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Figure 3. The radial distribution of X-ray surface brightness for a 
MOND isothermal sphere (T=2.5 keV) compared to observations 
of the Coma cluster (T=8.6 keV). The dashed curve is the /3 
model fit (Reiprich 2000) to the observations. 



is shown in Fig. 3 compared to that observed for the Coma 
cluster (Briel, Henry & Bohringer 1992). These observations 
are well fit by the traditional /3-model (Cavaliere & Pusco- 
Femiano 1976, Sarazin 1988), 



/(r) 



1+1- 



(6) 



shown by the dashed curve in Fig. 3. Here we see that the 
MOND isothermal gas sphere provides a poor representation 
of the actual surface brightness distribution. Moreover, the 

temperature of this best fitting MOND isothermal sphere is 
about 2.5 keV whereas the actual temperature of the Coma 
cluster is in excess of 8 keV. 



2.2 MOND /3-modeIs 

The radial dependence of electron number density which 
produces the X-ray intensity distribution described by eq. 6 
is 

-1.5/3 



Tie = Uo 



1 + 



(7) 



(Cavaliere & Fusco-Femiano 1976). Rather than assuming a 
constant temperature and solving eq. 1 for the density dis- 
tribution, one may alternatively take eq. 7 as the density 
distribution and solve for the radial dependence of the tem- 
perature. It is necessary to specify 13 (typically 0.6 to 0.7 for 
clusters), a central electron density, Uo (again taken to be 
0.006 cm^'^), a core radius, rc (ranging from 50 to 300 kpc 
for clusters), and a central gas temperature, To- Then eqs. 1 
and 2 may be solved for the run of temperature to an outer 
boundary, usually taken to be where the gas density falls to 
some fixed multiple of the mean cosmological density. For a 
given central temperature, the run of temperature is com- 
pletely determined by the core radius. There is one specific 



Figure 4. The projected emission-weighted temperature of a 
MOND 13 model as a function of projected radius (dashed curve). 
The radius is given in terms of Tmax where the gas density has 
fallen to 10~^* g/cm'^. This is the most nearly isothermal (3 
model. The solid curve shows the projected emission-weighted 
temperature for the most nearly isothermal two component model 
described in Section 3. 



value of Tc which minimizes the temperature gradients. For 
smaller core radii, the temperature rapidly increases toward 
the boundary (the models are very far from isothermal); for 
larger core radii, the temperature decreases to zero before 
the boundary is reached. Because clusters of galaxies are ob- 
served to be near isothermal, these MOND /3-models with 
the smallest temperature variations are taken to be appro- 
priate for clusters. 

The emission- weighted temperature of this most nearly- 
isothermal model as a function of projected radius is shown 
by the dashed curve in Fig. 4 for Uo — 0.006 cm"'', [3 = 
0.62, and To = 6.5 keV. The required core radius is 436 
kpc. It is evident from Fig. 4 that the model still deviates 
significantly from pure isothermal, with a roughly a factor 
of two variation around the central temperature. Moreover, 
compared to actual clusters, the implied core radius for this 
temperature is too large by a factor of two. 

For these near isothermal /3 models, 1 determined the 
mean emission weighted temperature, the total gas mass 
within the cut-off radius and the X-ray luminosity within 
the 0.1 to 2.4 keV band. The resulting gas mass-temperature 
relation is almost identical to that of MOND isothermal 
spheres (the dashed curve in Fig. 1) and thus again much 
larger than the observationally inferred gas mass of X-ray 
emitting clusters 

The luminosity-temperature relation for these MOND 
/? models is shown in Fig. 2 by the dotted curve. Unsurpris- 
ingly, this nearly coincides with the calculated relation for 
MOND isothermal spheres and is clearly an equally poor 
description of reality. The basic problem is that, for both 
MOND isothermal spheres and MOND f3 models, the core 
radii are too large. It is evident that some ingredient is miss- 
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Figure 5. The accumulated mass distributions in the Coma clus- 
ter. The dotted curve is the enclosed Newtonian mass as a func- 
tion of radius; the solid curve is the enclosed MOND mass; the 
dashed curve is the gas mass inferred from the X-ray observations 



ing from these models; that an additional mass component 
must be added to decrease the gas core radius at a given 
temperature. 



3 RESOLUTION OF THE PROBLEM: TWO 
COMPONENT MODELS 

For a cluster with an observed density and temperature 
distribution, eqs. 1 and 2 directly yield M(r), the interior 
dynamical mass as a function of radius. In the Newtonian 
regime this is given simply by 



Miv(r) 



kT 



G frup 



dln{p) _^ dln{T) 



dln(r) dln(r) 



(8) 



With MOND, taking n{x) to be given by eq. 2c, the dynam- 
ical mass is 

Mn 



Mm = 



^/l + {ao/af 



where a is the "observed" acceleration 



1 dp 
p dr 



kT 



frUpT 



dln{p) _^ dln{T) 



dln(r) dln{r) 



(9) 



(10) 



Obviously, from eq. 9 in the limit of large accelerations 
(a >> tto) the MOND dynamical mass is equivalent to the 
Newtonian dynamical mass. 

We may apply these relations to the Coma cluster which 
has a density distribution well-fit by a /3-model with /3 = 
0.71, rc = 276 kpc, Uo = .0036 (Reiprich 2001), an ob- 
served radial temperature profile (Arnaud et al. 2001) and 
average emission weighted temperature of 8.6 keV. The re- 
sults arc shown in Fig. 5 which is the cumulative Newtonian 
dynamical mass (dotted line), the MOND dynamical mass 
(solid line), and the gas mass (dashed line). Here it is ev- 
ident that while the Newtonian dynamical mass continues 



to increase at radius of 1 Mpc, the MOND mass has es- 
sentially converged. However, the total MOND mass is still 
a factor of four times larger than the mass in gas alone. 
This discrepancy cannot be accounted for by the stellar con- 
tent of the galaxies which, assuming a rnass-to-light ratio of 
seven, amounts only to about lO'^^ Mq within 1 Mpc (The 
& White 1988). This is, in fact, the discrepancy pointed out 
by The and White- a discrepancy which can be resolved by 
increasing ao by a factor of 3 or 4 over the value required 
for galaxy rotation curves, or by admitting the presence of 
non-luminous mass which MOND does not remove. 

The density of this non-luminous component is roughly 
constant and contained within two gas core radii. In other 
words the missing mass is essentially present in the inner re- 
gions of the cluster as implied in the work of Aguirre, Shaye, 
and Quataert (2001). This suggests that, with MOND, clus- 
ters might be described by a two component model: a gas 
component with a density distribution given by the /3-model 
and a dark central component of constant density and a ra- 
dius of about two times the core radius of the the gas dis- 
tribution. 

In Coma, the central surface density of the dark com- 
ponent is about Ed = 240 M^/pc^, which is comparable 
to the MOND surface density of a^/G ^ 700 Mq pc^. 
This is the characteristic central surface density of MOND 
self-gravitating isothermal systems (Milgrom 1984). There- 
fore, to determine scaling relations, I assume that the non- 
luminous component is a rigid sphere having a central sur- 
face surface density equal to that in the Coma cluster and 
radius twice that of the gas core radius, as in Coma. Then 
the constant density of the second component is 

P'i = ^ (11) 



2rc 



and the total dark mass is 



(12) 



Evidently, with this assumption, larger clusters have a more 
massive dark component, but because the gas mass scales 
as rc'^, the dark to gas mass decreases with increasing core 
radius or temperature. 

The procedure followed is identical to that of the pure 
MOND /3 models described above: I assume a gas distribu- 
tion described by eq. with /3 = 0.62, rio = 0.006. Then for a 
given central gas temperature, I determine the core radius of 
the model for which the temperature gradient is minimized. 
The new aspect is the second component which makes its 
presence felt by contributing to the total cumulative mass 
(M(r) in eq. 2) and hence the total gravitational force. This 
has the effect of decreasing the core radius at a given tem- 
perature compared to the single component MOND /3 mod- 
els. The emission-weighted temperature of as a function of 
projected radius is shown by the solid line in Fig. 4 again 
for a model with a central temperature of 6.5 keV. This is 
the most-nearly isothermal model, and we see that the tem- 
perature gradients are much smaller than in the the single 
component /3-model. The total variation about the central 
temperature is less than 40%. In other words, isothermal (3 
models require, in the context of MOND, this second central 
component with roughly constant density. 

The gas-mass-temperature relation of such models is 
shown by the solid line in Fig. 1 which is evidently consistent 
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with the observations. The luminosity-temperature relation 
for these two component cluster models is shown by the solid 
line in Fig. 2. These models provide a reasonable description 
of the observed relation. This is due to the fact that, in the 
low temperature systems, a relatively larger fraction of the 
mass is not in gaseous form. It is also not in the form of 
luminous material in galaxies as the implied mass-to-light 
ratios would be too large. Agreement of MOND with the 
cluster scaling relations is achieved at the expense of adding 
unseen matter. 



4 MODELING INDIVIDUAL CLUSTERS 

Individual clusters may be described by such two compo- 
nent models. Given the parameters of a /? model fit to the 
X-ray emission from an specific cluster (i.e., 13, rio, and 
Tc), and a characteristic temperature for the entire clus- 
ter, I determine, via eqs. 1 and 2, the central temperature, 
To, and the density (or surface density, related by eq. 11) 
of the dark component which yields the observed emission 
weighted temperature and the observed core radius, rc of 
the ,3 model. In all cases the dark component is assumed to 
extend to 2rc. Uniqueness is ensured by requiring that the 
temperature gradients be minimized- i.e., these are again 
the most nearly isothermal models. The fitting parameter is 
the density of the central dark component which yields the 
total dark mass via eq. 12 above. 

Table 1 lists the clusters which have been modeled in 
this way along with the the parameters of the /3-model 
fit and the mean temperature, all from the compilation of 
Reiprich (2001) with cluster properties scaled to h=0.7. The 
the required central surface density of the dark component, 
Ed, is given along with the total gas mass (out to the cut- 
off radius) and the total mass of the dark component. The 
enclosed Newtonian dynamical mass is also given. 

This sample was chosen to include a number of objects 
with significant inferred cooling fiows (such as Abel 1689 and 
2029). These clusters are characterized by relatively small 
core radii (rc < 100 kpc) and large central electron densities 
{rio > 0.01 crri^'^). The sample also includes clusters at the 
other extreme- large core radii (rc > 250 kpc) and low cen- 
tral electron densities {uo < 0.005 cm~^) with no inferred 
cooling flow (such as Abel 119 and 2256). 

In general, MOND reduces the classical Newtonian dis- 
crepancy for clusters of galaxies. For this sample, with 
MOND the ratio of the total dark mass to the gas mass 
is 1.60 ± 1.7. The Newtonian dark mass-to-gas mass ratio 
is 7.14 ± 2.7. While this is a significant reduction, it is also 
clear that MOND does not fully resolve the mass discrep- 
ancy in clusters. Moreover, it should be noted that the ratio 
of masses of the dark-to-hot gas components in the central 
two core radii (where the additional component is required) 
can be as large as 10, as was also pointed out by Aguirre, 
Schaye & Quataert (2001). This rigid component cannot be 
stars of a normal population because the mass-to-light ratio 
within two core radii would still be, on average, in excess of 
about 50. The principle question is whether or not this dark 
component is a fundamental problem or can be accommo- 
dated in the context of MOND. 

It could be that a high M/L population of low mass stars 
or sub-stellar objects is deposited in the central regions of 




log(dM/dt) (Mg/year) 



Figure 6. A log-log plot of the fitted dark mass (1O"M0) vs. the 
cooling mass inflow rate inferred from central gas densities and 
temperatures (White, Jones & Forman 1997). The crosses are 
those clusters with a small mass discrepancy {M^/Mg < 1.5), 
and the solid points are the objects with large discrepancies 
{Md/Mg > 1.5). 



clusters as a result of cooling flows. Fig. 6 is a plot of the total 
mass in the dark component vs. the cooling rate as estimated 
by White, Jones & Forman (1997). The solid points are those 
clusters with a large discrepancy: Md/Mg > 1.5. The crosses 
are clusters with a smaller discrepancy: Md/Mg < 1.5. We 
see that there is no obvious correlation between the total 
meiss of the dark component and the cooling rate— especially 
for those clusters with the largest discrepancy. On the other 
hand, in Fig. 7 we see a plot of the surface density of the dark 
component vs. the mass deposition rate. Here there does 
appear to be a correlation. This is because those clusters 
with the highest central dark matter densities are not the 
clusters with the largest mass discrepancies. 

It is unclear if this apparent correlation between surface 
density and mass deposition rate is significant. The mass 
deposition is not actually observed but calculated from the 
central gas density. Those clusters with large inferred cooling 
flows are clusters with high central gas dousitics and small 
core radii. But it precisely these clusters which require a 
large surface density of dark matter to produce the small 
core radius. While it may be the case that cooling flows 
contribute to the dark component of those clusters with the 
largest central density of dark matter, it is also evident from 
Fig. 6 that this cannot be the explanation for the discrep- 
ancy in clusters with the largest dark mass problem. These 
tend to be the clusters with low central gas densities and 
low inferred dark matter densities- but large core radii. In 
the next section, I consider another possibility: particle dark 
matter in the form of massive neutrinos. 
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Table 1. Two component model fits to observed clusters 
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42 
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0.99 
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5.2 
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5.0 
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4.06 
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All /3-model parameters are from the compilation of Reiprich 
(2001). In columns 8 and 10 the masses are those enclosed 

within a radius (typically 1 to 3 Mpc) where the gas den- 
sity has fallen to 250 times its mean cosmic value (assuming 
^5 = 0.04). In column 7 a blank entry means that there is no 
estimate for the cooling mass inflow rate. 



5 NEUTRINOS AS CLUSTER DARK MATTER 

The most well-motivated form of particle dark matter con- 
sists of neutrinos with finite mass. Aspects of the observed 
fluxes of atmospheric and solar neutrinos provide strong ev- 
idence for neutrino oscillations and hence non-zero neutrino 
masses (Gonzales- Garcia & Nir 2002). The fact that the 
number density of neutrinos produced in the early universe 
is comparable to that of photons then implies that there is 
a universal dark matter sea of neutrinos. The contribution 
to cosmological mass density would be il^h^ = Y] rrivf /94 
ev where the sum is over neutrino types. 



The neutrino oscillation experiments do not provide in- 
formation on the actual masses of neutrino species but on 
the square of the mass differences. These are small, such 
that, the largest mass difference, suggested by the atmo- 
spheric oscillations, is Am « 0.05 ev. If rriv ~ Am then the 
three active neutrino types would have no significant cos- 
mological mass density (« 10~^) and could not contribute 
to the mass budget of any bound astronomical object. But 
another possibility is that m„ >> Am and that the masses 
of all three active types are nearly equal. In this case an 
upper limit to the masses is provided by an experimental 
limit on the mass of the electron neutrino, i.e. 2.2 ev at 90% 
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Figure 7. The inferred cooling mass inflow rate in the sample 
clusters (as in Fig. 6) vs. the central surface density of the dark 
component in units of 100 Mq/pc'^. 

c.l. (Groom et al. 2000) If it were the case that the electron 

neutrino mass were near 2 ev, then neutrinos would consti- 
tute a significant fraction of the cosmic density {0,^, ~ 0.13 
for h=0.7). 

However, neutrinos of this mass could not contribute to 
the mass budget of galaxies. This follows from a classic argu- 
ment by Tromaino & Gunn (1979) based upon conservation 
of the phase space density of the neutrino fluid. Relativistic 
neutrinos are created with a maximum phase space density 
of (27r7i.)~^ per type including anti-neutrinos (this is a fac- 
tor of two less than the absolute limit implied by quantum 
mechanical degeneracy). In subsequent evolution of the neu- 
trino fluid involving gravitational instability and collapse, 
the final phase space density cannot exceed this value. This 
provides a relation between the final density of neutrino dark 
matter and the velocity dispersion of the system; with three 
types 

< (4.8 X lO-^^)(l^)^(Tfcev)^g/cm^ (13) 

Equivalently, for virialized systems, this may be written as 
a relation between the effective core radius of a dark halo 
and its velocity dispersion; this is, roughly, 

^ „ _/2ev\ VlOOOkm/s\ 5 ,^ „_/2ew\2,^ ,_i 

r-d > 0.5 '— Mpc w 0.7 (Tkev) * 

\ rriiy ) \ Or ) \ rriv J 

(14) 

Clearly for objects with the required velocity dispersion of 

galaxy halos (25km/s < (7^ < 200km/s), neutrinos in the 
mass range of one to two ev could not possible cluster on 
sub-Mpc scales. However, it obviously would be possible for 
such neutrinos to contribute to the mass budget of large 
clusters. 

Fig. 8 shows the the fitted matter density of the dark 
cores of the clusters of gala^Kies listed in Table 1 vs. the 
gas temperature of the clusters. The solid points are the 
clusters with significant mass discrepancies {Md/Mg > 1.5) 



Figure 8. A log-log plot of the fitted central density of the dark 
component in the sample clusters vs. the temperature. The solid 
line is the relation between the maximum neutrino density and 
the temperature imposed by phase space constraints for three 
neutrino types all with mass 2 ev (eq. 13) . The dashed line is the 
same when the neutrino mass is taJien to be 1 ev. The solid points 
are those clusters with a large discrepancy (M^/Mg > 1.5). 



and the open points the clusters with lower discrepancies 
{Md/Mg < 1.5). The solid and dashed lines show the rela- 
tion between maximum neutrino density and the tempera^ 
ture implied by eq. 13, for neutrinos of 2 ev and 1 ev respec- 
tively. For neutrinos of 2 ev, this relation appears to form an 
upper envelope for the clusters with the largest discrepan- 
cies. However, due to the extreme sensitivity of this density 
limit on neutrino mass (eq. 13), neutrinos with mass as low 
as 1 ev could not comprise the dark component. 

I have calculated the structure of a self-gravitating 
system supported by the degenerate pressure of neutrinos 
(a neutrino star). In the non-relativistic case there is a 
pressure-density relation of the form = KpJ"!^ where 
/sT = 5.5 X 10^^(m,./2eti)~3, and density is in g/cm^. In- 
serting this relation in eq. 1, I determine the run of density 
in a system with a given central velocity dispersion. Such ob- 
jects (n=1.5 polytropes) have nearly constant density out to 
a finite radius where the density drops to zero. As for white 
dwarfs, there is a mass-radius relation. This is shown by the 
solid curve (viv = 2 ev) and the dashed curve (jrii, = 1 
'^P^'ev) in Fig. 9 where I have assumed that the central velocity 
dispersion of the neutrinos is equal to the gas velocity disper- 
sion. The solid points show the mass and radius of the dark 
matter cores for those clusters with significant discrepancies 
{Md/Mg > 1.5) and the open points those with small dis- 
crepancies. The dotted line separates the MOND and New- 
tonian regime where we see that the mass-radius relation 
has different forms: r oc M^^^'^ in the MOND regime and 
r (X M~^^^ in the Newtonian regime. It is of interest that 
for Mv = 2 ev the range of radii and masses correspond in 
magnitude to the dark cores required in clusters. Moreover, 
we see that those clusters with large discrepancies lie gen- 
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Figure 9. The solid curve is radius-mass relation for self- 
gravitating objects supported by degenerate neutrino pressure, 
assuming three neutrino types of mass 2 ev. The dashed line is 
the same when the neutrino mass is taken to be 1 ev. The dotted 
line separates the Newtonian and MOND regimes and the points 
are relevant to the dark matter cores of the sample clusters. The 
closed points are those clusters with significant discrepancies as 
before. 



orally above the mass-radius relation as would be expected 
for neutrino clusters with density less than or equal to that 
given by eq. 13. Further, not even neutrinos of 2 ev appar- 
ently could contribute to the mass budget of those objects 
with low discrepancies. 

That neutrinos in the range of one to two ev make up 
the missing mass in clusters is a provocative possibility, and 
one which is within reach of experimental verification. The 
accelerator limits on the electron mass can soon be pushed 
to within a few tenths of an electron volt; if there is no 
positive detection, then active neutrinos cannot supply the 
dark mass of clusters. 



6 CONCLUSIONS 

What one may call, generically, "the dark matter problem" 
first became evident with radial velocity studies of rich clus- 
ters of galaxies (Zwicky 1933). The discrepancy between the 
visible and Newtonian dynamical mass, quantified then in 
terms of mass-to-light ratio, was more than a factor of 100. 
With the advent of X-ray observatories and the detection 
of hot gas in clusters, this discrepancy between dynamical 
and detectable mass was reduced to a factor of 10. Modified 
Newtonian dynamics reduces the discrepancy further, to a 
factor of 2 to 3, but it is clear that a discrepancy remains 
which cannot be explained by detected gaseous or luminous 
mass. It is also apparent that while the missing mass is pri- 
marily in the inner regions of clusters, it does extend beyond 
the core radius as defined by the gas density distribution. 

How serious is this remaining mass discrepancy for 
MOND? Formally speaking, it does not constitute a falsi- 



fication. If the dynamical mass predicted by MOND were 
generally less than the detected mass in stars and gas, then 
it would be a definite falsification, but this is not the case. 
More mass can always be found (as in the case of the hot 
gas), but it is difficult to make observed mass disappear. 

One might reasonably argue that the implied presence 
of undetected mass runs against the spirit of MOND, which, 
in the view of most people, was suggested primarily as a 
replacement for dark matter. But more generally, MOND 
should not be viewed simply as an alternative to dark mat- 
ter; the systematic appearance of the mass discrepancy in 
astronomical systems with low internal accelerations is an 
indication that Newtonian dynamics or gravity may break 
down in this limit. MOND primarily addresses this issue: 
is physics in the low-acceleration regime Newtonian? The 
success of MOND in explaining the scaling properties and 
observed rotation curves of galaxies suggests that it may not 
be. MOND does not rest upon the principle that there is no 
undetected or dark matter. Indeed, comparing the density 
of luminous matter to the baryonic content of the universe 
implied by considerations of primordial nucleosynthesis, one 
can only conclude that there is, as yet, undetected baxyonic 
matter, probably in the form of diffuse gas in the intergalac- 
tic medium. Moreover, it is virtually certain that particle 
dark matter exists in the form of neutrinos; only its contri- 
bution to the total mass density of the Universe is unclear. 

MOND would be incompatible with the wide-spread 
existence of dark matter which clusters on the scale of 
galaxies- cold dark matter. But MOND is not inconsistent 
with hot dark matter such as 2 ev neutrinos, which can 
only aggregate on the scale of clusters of galaxies- indeed, 
I have presented evidence that this may be the case. Neu- 
trinos, as particle dark matter candidates, are unquestion- 
ably well-motivated, both from a theoretical point-of-view 
(they definitely exist) and from an experimental point-of- 
view (they have mass). No conjectured CDM particle shares 
these advantages. While I do not wish to state that the dark 
matter in clusters is definitely in the form of 2 ev neutrinos 
(there is more than enough remaining baryonic matter to 
make up the missing mass), there are indications that point 
this way. The largest discrepancies are found in the clus- 
ters with the largest core radii as would be the case with 
neutrinos. The indicated densities of dark matter are com- 
parable to the maximum possible density of 2 ev neutrinos. 
The radius-mass relationship for self-gravitating degenerate 
neutrino objects forms an envelope for those clusters with 
large discrepancies. 

The fact remains that there exists an algorithm, 
MOND, which allows galaxy rotation curves to be predicted 
in detail from the observed distribution of matter, and it is 
for these systems that the kinematic observations are most 
precise. This fact challenges the current CDM paradigm, and 
demands explanation if dark matter lies behind the discrep- 
ancy. The factor two remaining discrepancy in clusters is 
less challenging for MOND, particularly given that MOND 
makes no claims about the full material content of the Uni- 
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